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Abstract: 

Radiative corrections to both perturbative and non-perturbative contributions are added 
to existing calculations of the Isgur-Wise function ^jw- To this end, we develop a method 
for calculating two-loop integrals in the heavy quark effective theory involving two dif- 
ferent scales. The inclusion of 0{as) terms causes to decrease as compared to the 
lowest order result and shows the importance of quantum effects. The slope parameter 
violates the bound given by de Rafael and Taron. 
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1 Introduction 



The heavy quark effective theory (HQET) |T|, |2[ now seems estabhshed as a useful tool 
for investigating the low energy regime of QCD in the limit of infinitely heavy quarks. 
New flavour symmetries appearing in the Lagrangian allow for a considerable simpli- 
fication of the treatment of transition amplitudes of heavy particles. Since the heavy 
quark decouples from the light degrees of freedom, the (heavy) meson to (heavy) meson 
amplitudes, for instance, are completely determined by a single universal function, the 
Isgur-Wise function ^/vf 0. ^/^y describes the dynamics of the light quark cloud and de- 
pends on the "momentum-transfer" y := v-v', where v and v' are the four- velocities of the 
two heavy mesons, respectively. It satisfies the normalization condition C,iw{.y)\y=i ~ ^ 
and one can show that the slope at the normalization point y = 1 is determined to be 
p2 := - d^iw{y)/dy\y=i > 1/4 0. In addition, it has been claimed that < o.45 [g. 
As far as we know, the latter bound is in contradiction to almost all parametrizations to 
be found in the literature P . In view of this discrepancy it seems worthwile to address 
the problem by means of HQET sum rules which were first developed for the (heavy) 
meson-to-vacuum transition amplitude 0, §|, and yield nice agreement with lattice cal- 
culations. Since the former calculations stress the importance of radiative corrections, it 
is of great interest to include them in the Isgur-Wise function as well, whereas they have 
been omitted from existing sum rule calculations of [H ii Hof. In this letter we now 
present a calculation of the Isgur-Wise function in the framework of HQET sum rules 
including radiative corrections up to 0{as). 

Our paper is organized as follows: in Sect. 2 we develop a method for calculating 
two-loop integrals in the HQET involving two different scales and velocities as a series in 
r := {y — 1) / {y + 1) and give our result to 0{as) for the three-point function relevant to 
the determination of dw- In Sect. |^ we derive a sum rule for the renormalization group 
invariant quantity Finally, Sect, ^is devoted to the discussion of our results. 

2 Two-loop integrals in the HQET 

In order to obtain the Isgur-Wise function within the framework of HQET sum rules we 
first need to evaluate the three-point function 

f{uj = q- V, u' = q'- v', y) (v + v')^ = t'^Jd'^x d% e^'^'^-'^'-^ { | T Mx)V,iO) JUy) | ) , (1) 

with the HQET-currents J5 = qi'j^h^ and = h^'J^ih^', where h^j denotes a heavy quark 
field in the HQET with four-velocity v and q denotes a light quark field. 

In this section we shall be mainly concerned with the 0{as) (two-loop) correction, 
f pelt; to the perturbative part of the three-point function. Throughout this section we 
will drop the subscript pert, for no confusion can arise. After some algebra, one can write 



- 1 - 



combination of two-loop integrals of the form 

1 \'^/ 1 \"/ 1 / 1 X"^/ 1 



Vfc^/ \{l — ky J \v ■ I + uj J \v' ■ I + uj' J \v ■ k + uj J \v' ■ k + uj' 

(2) 

where D = 4 — 2e is the space-time dimension. We have not succeeded in writing the 
result of these integrations in a closed form convenient enough for numerical analysis. 
We have overcome this difficulty by expanding v' and v around {v' + v)/2, thus obtaining 
an expansion of the integrals in powers of r = (?/ — !)/(?/ + 1). For the numerical analysis 
of Sect. ^, it will be sufficient to truncate the series at order r^. Since this expansion 
technique can certainly be useful in other 0{as) calculations, we shall present it here 
briefly. 

We begin by introducing the variables f+ = (f + v')/2 and = {v — v')/2. Because 
of f ^ = f'^ = 1, we have f + ■ f _ = 0. We now expand all v- and f '-dependent factors in 
eq. (ID in powers of f _ ■ A; or t>_ ■ /. To get the required accuracy, we need to keep only 
the first five terms of the expansion. We end up with integrals of the form 

r d^ld^k{v_-iy{v_-k)'3 
J (/2)«(A;2)&((/ - kyy{v+ ■ I + cu)"(i;+ ■ / + uj'y{v+ ■ k + a;)™(v+ ■ k + uj'Y ' 

This expansion removes the main source of trouble, namely the presence of two different 
velocities in the denominators. We next scale the loop variables as 



r ^ ^vl , F ^ F ^vl, (4) 
and introduce the notation v^^ = vl/^Jv^, so that -O^ = 1. Eq. (H) now reads 

■I _^ y\{2a+2b+2c-2D-p-q)/2 



d^l d^k ■ ■ • Ifxpky^ ■ ■ ■ k^^-^ 



X 



where the indices in parenthesis are to be symmetrized. The integral in the last equation 
must be of the form Y.k-^kiXklfi^- uq where (7^)^^...,^^ are totally symmetric tensors built 
up from g^y and Note that those tensors containing -0^, can be ignored, since 
v-V-=Q. As a consequence, J(p+q=odd) = 0, whereas one can show that J{p+q=2n) = 0{r'^). 
Thus, only the cases p + q = 0, p + q = 2 and p + q = 4 need be considered. The 
coefficients Ak can be obtained by the usual procedure of contracting all indices of the 
integral in eq. with the appropriate number of metric tensors, g'^", and velocities, v". 
This will provide a set of equations for Ak which can be solved in terms of scalar integrals 
of the form 



I^'\a,b,c,p,q)= Jd^'ld^ki^—^ i^— 



{k-lYj \v-k + ujj U-/ + tu(') 



(6) 
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Those integrals I depending on only one scale uo have been evaluated in Only 
I'{a,b,c,p,q) deserves some comments. For p > 1 and/or g > 1, I'{a,b,c,p,q) can be 
obtained from I'{a,b,c, 1,1) by simply taking derivatives with respect to u and/or uj'. 
To evaluate /'(a, b, c, 1, 1) one can make use of the following master equation: 



c-aH (2a + c + p- D) 

UJ 



I'{a,b, c,p,q) 



(7) 



c(l-ii)C+(A- 

UJ 



B-) + aA+{B- - C-) + —pP^Q- 

UJ 



I'{a,b, c,p,q), 



where the action of is given by A^I'{a,b,c,p,q) = I' {a ± l,b,c,p,q), with similar 
definitions of 5^, C^, and Q^. We have deduced eq. (|^ by following the "integration 
by parts" procedure also used in [|TT]. By means of eq. (0), one can reduce any 



/'(a, 6, c, p, g) appearing in the calculation of f^^-* to cases where at least one of the 
arguments a, b, c, p, or q vanishes. The cases p = 0, q = 0, c = can be evaluated 
using the formula given in For the cases a = and 6 = we make the change 



k,l 



{I + k), and use the results of [TTI] as well as 



-1 
1^ 



UJ 



UJ 



UJ + V ■ kj \uj' + V ■ k 



it: 



D/2 



T{2a + p + q- D)T{D/2 
r(a)r(p + g) 



X 



UJ 



UJ' 



'-2uj 



l\D-2a 



2Fi{p + q + 2a- D,p,p + q,l 



UJ 



UJ' 



where 2Fi{ci,b,c,t) is the hypergeometric function. 

Using a similar technique, it is straightforward to obtain the lowest order contribution, 
T^^K Up to 0{as), the full three-point function can be written as 



~ , div fin t-"- e ^ o div II , div I fin 



;(o) 



(1) 



(9) 



where the dependence in e is given explicitely. It is important to remember that the 
0(ea°)-terms in eq. (U), i.e cannot be neglected since they give a finite contribution 
of order to F through the external current renormalization. The terms F^*?!, F^j^ jj are 
polynomials in tu, uj' and r whereas the other contributions in eq. contain logarithms 
and dilogarithms as well as negative powers of uj — uj'. By expanding these functions in 
powers of uj' — uj, one can check that the limit cu' — > of F exists. Actually, F is analytic 
in uj' — UJ for uj, uj' < 0. Note that the UV-divergent portion of f i.e. f j has got 
a non-polynomial structure in uj, uj'. As in ordinary QCD, these non-localities disappear 
upon renormalization of the currents J5 and V^. The corresponding renormalization 
constants, Zhl and Zhh have been computed by 0, Q to be 



1 + 



as 
47re 



1 + 



as 
47re 



32 64 2 3, 



(10) 
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We have checked that the UV-divergent part of the three-point function of the renor- 
mahzed currents is a polynomial in the external variables a;, uo' , as it should. This is a 
nontrivial check of our result. Hereafter, F will stand for the renormalized three-point 
function. 

As usual in the sum rule approach, one would like to express F as a dispersion integral, 

i.e. 



1 

vr^ Jo 



ds [ ds'- — ' ' — - + subtraction terms 
Jo [s — uJ){s' — to') 



'111 



with psi^u, u', r) 



Im <^ Im F(u;, s' + iO, r] 



(12) 



Applying eq. ([12|) to our result for f we obtain 

1 d^ 
- —p3{s, s',r)=Qo{s, s', r)—6{s - s') 



(13) 



d'^ d'^ 
+Qi{s, s', r)-j^log \s-s'\+Q2{s, s', sgn(s - s') log ], 

where Qo, Qi, Q2 are polynomials in s, s' and r. By substituting the explicit result sum- 
marized in eq. (0) in the partial derivative of eq. (pH]) with respect to 00 and 00', we have 
checked that one recovers the full expression for d'^r{uj,uj',r)/dujduj' computed directly 
by differentiating our original result for r{uj,uj',r). The differentiation is necessary in 
order to eliminate the subtraction terms in eq. (|TT|). 

Using the explicit form of eq. ( [TB| ) it is quite straightforward to compute a quantity 
that plays a central role in the next section: the Borel transform of F, B^'^^B'^'^'^T. One 
just has to evaluate the integral 



BM.B'M.f = _ 1 r ds r ci.'p3(^,^',2/)e-(^+^')/^^ 

TT^Ml Jo Jo 



(14) 



where M3 is the Borel parameter for the three-point function. The final result is a very 
simple expression which can be cast into the following form: 



ttWI Jo 1 8 



17 31og(2s//i) 



+ C(2) 



f 3 fag 



355 41og2 13, , , , 

+ ^^ + ^log(2V/^)-2C(2) 



72 



6 



+ 



11873 32, ^ 109, ,^ , , 

-— log2-— log(2s//i) + C(2) 



3600 15 



60 



(15) 
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3 Evaluation of the Isgur-Wise function 



We are now in a position to derive a sum rule for the Isgur-Wise function including the 
lowest order radiative corrections. Sum rules for the meson-to-vacuum amplitudes in the 
HQET were introduced by |] and applied to the Isgur-Wise function by P, ||, . We 
will not go into details about the method of HQET sum rules, for which we refer to |Q , 
but simply review the main points. 

One can saturate eq. (|lj) with physical states of the HQET yielding 

had I T^cont 

_ ( I Js I P'{v') ) ( P'{v') I V, I P{v) ) ( P{v) I J5 I ) 
4(Am - uj){Am - 00') 

+ contributions of higher states, (16) 



r{v + v'), = rj'^'^ + r^ 



where Am = tjim — itiq ^ 0.5 GeV is the difference between the meson and the quark 
mass in the heavy quark limit (HQL). The matrix elements in the numerator can be 
expressed in terms of the leptonic decay constant in the HQL, /, 

{0\.l\P) = f, (17) 

and the universal Isgur-Wise function 

{P'\V,\P)=^jwiy)iv + v'),. (18) 

On the other hand, F can be calculated by means of the operator product expansion 



(OPE) a la |15|, away from the physical region, i.e. for uj, uj' <^ 0. Both expressions are 



matched after application of the Borel transformation defined in (|1J) which suppresses 
the contributions of both higher operators in the OPE and higher states in ([T6|) . 

Due to the scale-dependence of B^^ j^'^s -^g shall derive a sum rule for the renor- 
malization group independent (to two-loop accuracy) quantity 

TT 



f ^ §Ah B'Kh f (/i)a,(/x)-(27o^"+^o"")/(2/3o) [ 1 _ ^iM |2A^^ + A^^} ] (19) 



with A = f 21 _ I) and 7"-(-^) = (70""^"'^ + f^7r^"'^ + ...)> (20) 
where 7^-'^ {l^^) denotes the anomalous dimension of the effective current with two 



heavy ( one heavy and one light) quarks ||TT], |T3|, |Tj]. [3q and [3i are the lowest order 
coefficients of the usual /5-function of QCD. Following the analysis of [|], we choose 
/i = M3 which is the natural scale for the Borel transformed expression. 

As for the OPE, we take into account operators up to dimension 5 or, to be specific, 
the quark condensate ( qq ) and the mixed condensate ( qagGq ) . The gluon condensate 
gives a small contribution and only at 0(|/ — 1) 0. Likewise, the contribution of the 
four-quark condensate is expected to be small. Thus we write 

Qhh §'Ah f = fpert ^ f{3)Q^ ^ f^^^O^ + ... (21) 
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with the renormahzation group invariant quark and mixed condensate Os and O5, re- 
spectively, as defined in f^. The perturbative contribution follows from (|T5|): 



TT 



17 , 



2 4 i M3 9 



16 , 49 
— In2 

9 54 

^+0.3 



8 , 197 

— In2 + 

15 600 



-0.04 



(22) 



Here we have factored out an overall factor l/(?/ + 1)^ which appears in the exact expres- 
sion for the bare loop diagram. Furthermore, we have summed up the coefficients in front 
of the logarithm to get the anomalous dimension as specified by the scale-dependence 
of r. The numbers under the underbraces show the good convergence of the series in 

The lowest order expression for F^^^ is simply given by — 1/(4M|), whereas for the 
0{as) corrections we take the expression given by |p. Radiative corrections to F^^-* have 
not been calculated so far, so we use the lowest order expression F^^-* = (2?/-|-l)/(192M|). 

Using the usual argument of quark-hadron duality, the contribution of higher states to 
(16) is modelled by the perturbative contribution to the OPE (^) above certain threshold 



As, the so-called continuum threshold. Following the discussion of ||16[, we choose the 
threshold in such a way that undesirable contributions of P-waves or even higher states 
do not mix in. It was already noticed by [^, |[ that the slope of the Isgur-Wise function 
is quite sensitive to the choice of the continuum model and introduces a rather large 
uncertainty. This problem gets resolved by choosing the continuum model according to 
161 as 



^cont 



ds X Integrand of (22). 



(23) 



AS3 



Taking all together we get the following sum rule for the renormahzation group invariant 
Isgur-Wise function: 



iiwiy) 



4M|(5^^3^/A/3f(^) _f 



cont\ 



(24) 



where B^^^ SR denotes the corresponding sum rule for the square of the leptonic decay 
constant, and can be found in |p. Note that the exponential exp(— Am/Ms), that 
formally appears on the left hand side of (^) due to the Borel transformation, cancels 
against the corresponding factor in ]3^''^ SR and, hence as calculated in (p^ becomes 
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independent of Am. The sum rule (E^) automatically fulfills the normalization condition 



^iwiy) = 1 for any value of continuum threshold and the Borel parameter provided 

y=l 

we take the same value for As2 = Ass =: As in both the numerator and denominator 
and take M3 = 2M2 =: 2M. This equivalence will be used throughout the next section. 



4 Results and Discussion 

We now turn to the numerical analysis of the sum rule eq. (^41) . The numerical values 
of the condensates we use are {qq) = (—0.24 GeV)'^ and {qcgGq) = (0.8 GeV^)(gg) 
at the scale of 1 GeV. In Fig. |I| we show C,iw for different values of As and M as a 
function of y. For phenomenological applications, such as the decay B D*ev, it is 
sufficient to know the Isgur-Wise function within the range 1 < y < 2. In this region, the 
sensitivity of eq. (0) to the continuum threshold and the Borel parameter is found to 
be quite small and to amount to at most 6% for y = 2. In the evaluation we have taken 
1.1 GeV < As < 1.4 GeV and 0.5 GeV < M < 1.0 GeV as suggested by the analysis 
of the sum rule for / |^. The resulting curves may be parametrized by a second order 
polynomial in {y — 1) as 

iiwiy) = 1 - (0.54 ± 0.01)(i/ - 1) + (0.17 ± 0.01)(y - 1)2 (25) 

where the errors reflect the uncertainty due to M and As. The result is mainly determined 
by the perturbative term which contributes between 50% and 85% at y = 1, depending 
strongly on M, and increasing with y. 

In Fig. we depict the effect of neglecting radiative corrections which is less than 
10% for y < 2. It is, however, clearly seen that the inclusion of radiative corrections 
lowers the values of ^iw- This is not unexpected an effect, since the coefficients of powers 
of {y — 1) in eq. (0) are quite small and overcome by the ?/-dependent terms in front 
of the logarithm. On the other hand, the large constant term ~ as/n {17/3 + 47c'^/9), 
that was attributed to Coulombic corrections in [§], is cancelled to a large extent by the 
corresponding term in i?^^^ SR. In other words: the radiative corrections are mainly 
determined by the one-gluon exchange between the heavy-quark lines, i.e. radiative cor- 
rections to the weak vertex, and thus pure quantum effects (and not classical ones such 
as Coulombic corrections). 

In order to determine the slope parameter of the Isgur-Wise function, = — 
{dC,iw{y)) / idy)\y^^, "we scale C,iw to a physically meaningful scale as 

^iw{y, m) = [a.(m)]-o-"/(2*) ^ oJ^^hh^j ^^^^y^ (36) 

To be specific, we choose fh = mBTnD/ {ttlb + m/j) ^ 1.4 GeV, the harmonic mean of the 
masses of the B- and D-meson. 
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From the above formula we get = (0.84 ± 0.02) where the error stems from the 
uncertainty in the choice of parameters. This value is nearly twice as big as the upper 
bound Pmax = 0-45 given by @]. Disregarding all radiative corrections but the overall 
leading-log scaling factors, we get = (0.78 ± 0.02) which still is far away from p^ax- 
We see no possibility to get such a small value for the slope parameter, but rather agree 
with the existing calculations and fits which yield p^ ~ 1.0 — 1.5. The resolution of 
this disagreement surely requires further investigations. 
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Figure Captions 



Figure 1: ^iw as a function of?/, eq. (H, for 0.5 GeV < M < 1.1 GeV, 1.1 GeV < As < 
1.4 GeV. The spread of the hues reflects the uncertainty due to the choice of parameters. 



Figure 2: Sohd hue: ^iw as a function of y, eq. (|2|), at M = 0.9 GeV, As = 1.23 GeV. 
Dashed hne: without radiative corrections. The influence of 0{as) corrections causes 
a stronger falling-off of the form factor mainly due to the corrections to the weak vertex. 
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